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Quantum simulation of a Hamiltonian H requires unitary operator decomposition (UOD) of its 
evolution operator, (U — exp(-iHt)) in terms of experimentally preferable unitaries. Here, using 
Genetic Algorithm optimization, we numerically evaluate the most generic UOD for the Hamiltonian, 
DM interaction in the presence of Heisenberg XY interaction, Hdh- Using these decompositions, 
we studied the entanglement dynamics of Bell state in the Hamiltonian Hdh and verified the 
entanglement preservation procedure by Hou et al. [Annals of Physics, 327 292 (2012)]. 



I. INTRODUCTION 

Algorithms with exponential speedups over classical 
counterparts [TJ [2] , simulation of quantum systems [3H5] 
and testing basic principles of quantum mechanics [6j [7] 
makes quantum Information processing (QIP) and quan- 
tum computation intensively investigated fields of physics 
over last decade. The idea of simulating quantum sys- 
tems in a quantum computer was proposed by Feynman 
[8] in 1982 and is one of the most important practical 
application of the quantum computer. Quantum simu- 
lation has the potential to revolutionize the physics and 
chemistry and draws attention recently by solving prob- 
lems like - molecular Hydrogen simulation [4 j, calcula- 
tions of thermal rate constants of chemical reactions [5] 
and quantum chemical dynamics [9]. 

Dzyaloshinsky-Moriya (DM) interaction is an 
anisotropic antisymmetric exchange interaction arising 
from spin-orbit coupling [TO] [TT]. It was proposed by 
Dzyaloshinski to explain the weak ferromagnetism of 
antiferromagnetic crystals (a- F 62 O3, MnCOi) [TO]. 
DM interaction is crucial in the description of many 
antiferromagnetic systems [T2HT4] and is important in 
the entanglement properties of the system. Here we 
present a generic unitary operator decomposition which 
will help to simulate the Hamiltonian - DM interaction 
in the presence of Heisenberg XY interaction - in a two 
qubit system with almost any basic interaction between 
them. 

Long-lasting coherence and high fidelity controls in nu- 
clear magnetic resonance (NMR) are ideal for quantum 
information processing. Experimental implementation 
quantum algorithms (Deutsch-Jozsa algorithm, Grover's 
search algorithm and Shor's algorithm of factorization), 
testing basic principles of quantum mechanics (nohiding 
theorem [6] and Leggett-Garg inequality [7]) and quan- 
tum simulation (hydrogen molecule [4 and system with 
competing two and three Body interactions [15]) were 
performed in Liquid state NMR (LSNMR). 

Genetic algorithms (GA) are stochastic global search 
method based on the mechanics of natural biological evo- 
lution [16]. It was first proposed by John Holland in 
1975 fT7 . GA operates on a population of solutions of 
a specific problem by encoding the solutions to a sim- 
ple chromosome-like data structure, and applies recom- 



bination operators. GAs are attractive in engineering 
design and applications because they are easy to use and 
are likely to find the globally best design or solution, 
which is superior to other design or solution [18]. Here 
we used Genetic algorithm optimization for solving UOD 
for generic DM Hamiltonian with Heisenberg-XY inter- 
action. 

Section H] deals with theoretical discussion of DM 
Hamiltonian simulation followed by experimental imple- 
mentation in Section Hill 

II. THEORY 

DM interaction in the presence of Heisenberg XY in- 
teraction H(J, D) is, 

H(J,D) = J( y CFi x CF2 x J rCFiyCF2y) J rD(cFi x CF2y-CFiyCF2x), (1) 

where J and D respectively represents the strength of 
Heisenberg and DM interactions. 

Experimental simulation of H ( J, D) (Eqn. |lj in a 
quantum system (with Hamiltonian H sys ) requires UOD 
of evolution operator /7(J, £>,£), 

U(D, J, t) = exp(-iH(J, D) x t), (2) 

in terms of Single Qubit Rotations R n (6,(f)) (6 angle ro- 
tation along cj) axis on n th spin), 

R n (0, (j)) = exp(-i0/2 x [Coscj) cr nx + Sin<j> cr ny \), (3) 

and evolution under system Hamiltonian U sys (Eqn. |4|, 

U sys (t) = exp(-iH sys x t). (4) 

Without losing generality, Eqn. [2] can be written as, 

U (7, r) = exp(-i[(a lx a 2x + criyV2y) + 

l(cf lx CF2y - Cr ly (J2x)] t), 

where 7 represents the relative ratio of interaction 
strengths (7 = D / J) and r = J x t. 

Eqn. [5] forms the complete unitary operator for the 
Hamiltonian (Eqn. [I]) with 7 and r varies from to 00. 
We performed UOD for Eqn. [5] using Genetic algorithm 
optimization [19 j. In an operator optimization (as shown 
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FIG. 1. Fidelity profile of UOD given in Eqn. [7| 



in [19]), optimization is performed for a constant uni- 
tary matrix- corresponds to a single fidelity point. Here 
optimization has to be performed for a two dimensional 
fidelity profile generated by 7 and r. We name it as Fi- 
delity Profile Optimization (FPO). FPO for the present 
case is explained in following steps. 

In the first step, we performed Fidelity Profile Opti- 
mization with following assumptions - (a). the range of r 
is from to 15, (b). the range of 7 is from to 1 and 
(c). the system Hamiltonian (H sys ) is given by Eqn. [6] 



H sys — H zz — Jzz{?\z (J 2z)- 



(6) 



where J zz is the strength of zz-interaction. 

The optimization procedure using Genetic algorithm 
is explained in the Supporting information. 

The optimized UOD (Eqn. [7]) has seven SQRs and two 
system Hamiltonian evolutions. 



C/( 7 , r) = R 1 ( f , - f )R 1 ( f ,9 2 )R 2 (tt, tt) U zz ( f ) 

R\e 1 ,e 2 + i)R 2 {'K-e u Q)u zz {i) 

i? 1 (f^ 2 + 7r)i? 2 (f,f), 



(7) 



where 0\ and 62 (Eqn. [8| impart 7 and r dependence to 
UOD and U zz is given by Eqn. [4] with J zz x t = 7r/4. 



0i = [0.8423 - 0.3455 Cos(1.117 7)+ 
0.01806 5m(1.1177)] r 
2 = 1.345 exp(-0.873l7) + 1.796. 



(8) 



The fidelity [19] profile of UOD is shown in Fig. [T] The 
minimum point in fidelity profile is greater than 99.99 %. 
It should be noted that, the total length of UOD (Eqn. 
[7| is invariant under 7 and r (with the assumption - all 
the SQRs are instantaneous). 

For generalizing the assumption on r, we solved Eqn. 
9] numerically and find the period P(7) of U(7,r) (Eqn. 
10b. 



#(7, r + n x P( 7 )) = if (7, r). 



(9) 



P(7) = 3.008 7 3 - 6.627 7 2 - 0.1498 7 + 12.59. (10) 



Eqn. 10 has a maximum value of 12.59 at 7=0. Since the 
maximum value of period is less than 15 (FOP performed 
till r=15), UOD (Eqn. can be used for any value of 
r. Same argument can be used for extending the range 
of r to —00. 

In order to incorporate the range of 7 from to 00, we 
performed FPO for the operator Eqn. 11), 



U'(i, T ') = exp(-i[j f (a lx a 2x + cr l3/ a 23/ ) + 

{Vlx°2y — Vly&2x)\ T*), 

and the optimized unitary decomposition are, 

[/'( 7 V) = ^(1 , |)# 2 (| ,e 3 )Uzz{l)R\o 2 + 1 

R 2 (0i,0 4 )U zz (l)R 1 ^,^)R 2 r 



(11) 



5,0) 
2 » 6*3), 

(12) 



where 6\ • • • 64 (Eqn. [13| ) impart 7 and r dependence to 
UOD. 

= [0.09812 exp(-2.42 7 ) + 0.4023 exp(0.5524 7 )]r, 
0! = -61 + 3.142, 

6 2 = e- [1.242 exp(-0.96177) + 0.3546 exp(-0.11457)], 

3 = 1.259 exp(-0.9577) + 3.479 exp(-0.00877), 

4 = 1.256 exp(-0.9597) + 1.912 exp(-0.01667), 

(13) 

where 7' varies from to 1 and r from to 15. 

Eqn. [TT| sa tisfy the same periodicity relation as shown 
in Eqn. |10] and hence can use the same reasoning for 
extending r range from to +00. 

For 7 > 1, Eqn. [5] can be written as, 

U(n",r") = exp(-i[j" (ai x a 2 x + 0i y (T2y) + 



{?\x°2y - Vly°2x)] 



(14) 



where 7" = I/7 and t" = 7 x r. 



Eqn. [11] and Eqn. [14] are equivalent and hence UOD 
for Eqn. [5]can be shown as, 



C/( 7) r) = 



Eqn. 
Eqn. 
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if 7 < 1 
if 7 > 1 



(15) 



H. 



The UOD optimization given Eqn. [7] is based on 
H zz (Eqn. [6|). It can be generalized to almost 



sys 



any interaction by term isolation procedure by Bremner 
et al. [201 HI]. 

As an example consider the case, 



J(0- X CF X + CFyCFy + G Z G Z ). 



(16) 



The terms can be isolated from Eqn. 16 and is 



shown in Eqn. 17 



exp(-iJ zz a z a z t) = R 1 (it , z)exp(-iH sys t)R 1 (it , z) 



(17) 



exp(-iH sys t), 

where R 1 (tt^z) represents a 7T-SQR on spin 1 along Z 



axis. 



Combinig all the steps above forms most generic UOD 
of the Hamiltonian - DM interaction in the presence of 
Heisenberg XY interaction. 
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FIG. 2. C labeled Chloroform used for quantum simulation. 
13 C and 1 H act as qubits with zz interaction (Jc# =2 15.1 Hz) 
between them. 



III. EXPERIMENTAL QUANTUM 
SIMULATION 



We performed Quantum simulation experiments in a 
two qubit NMR system 13 CHCls (dissolved in Acetone- 
D6) (Fig. 2| with 13 C and 1 H spins act as two qubit 
system with scalar coupling (zz interaction- Eqn. [6| be- 
tween them. The system Hamiltonian is zz interaction 
(Eqn. [6J. We performed all the experiments in Bruker 
AV-500 spectrometer. 

Quantum computation experiments in NMR starts 
with (i). preparation of pseudo pure states [22-24 , (ii). 
processing the state by evolving under different average 
Hamiltonians [25] and (iii). read-out by quantum state 
tomography [26]. Here we studied the entanglement dy- 
namics (quantified by concurrence [27]) of a Bell state 
(Eqn. (181) in the Hamiltonian given in Eqn. [TJ 
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FIG. 3. (a). Entanglement (concurrence) dynamics of 13 C — H 
system under the Hamiltonian Eqn. Jll (b). Entanglement 
preservation experiment using Eqn. |20| Starting from sin- 
glet state, the concurrence sustains at 1 with the preservation 
procedure. 



Using the unitary operator decompositions shown in 
we have simulated the Hamiltonian if (7, r) 
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Eqn. 

for 7={0.33, 0.66, 0.99} and studied the entanglement 
dynamics of the singlet Bell state (Eqn. [l8| under these 
Hamiltonians. 

In experimental implementation i? n (#,(/>) is imple- 
mented by hard pulse [28] with suitable length (deter- 
mined by 0) along the axis <ft on 1 H or 13 C spin and 
U zz (9') is implemented by system Hamiltonian evolution 
for a time determined by 0' [25]. The experimental sim- 
ulation results (Fig. 3(a)) shows a good agreement with 
the theoretical simulation. Average experimental devia- 
tion (AED) in concurrence - calculated using the formula 
Eqn. [l9]- is 3.83%, which is mainly due to decoherence 
effects, static and rf inhomogeneities. 



1. Entanglement Preservation 

Hou et al. [29] demonstrated a mechanism for entan- 
glement preservation of a quantum state in a Hamiltonian 
of the type given in Eqn. [2] 

Preservation of initial entanglement of a quantum state 
is performed by free evolution interrupted with a certain 
operator O, which makes the state to go back to its initial 
state. The operator sequence for preservation is given in 
Eqn. [201 



OC/(7,r)OC/( 7 ,r)= J, 



(20) 



AED = ^ c^) 



(19) 



1=1 



where C es (i) and C ts {i) are the concurrence in exper- 
imental and theoretical simulations and n is the number 
of experimental points (here we performed simulation for 
n = 16). 



where O = 1 1 (8) (J2z • 

We performed entanglement preservation experiment 
for singlet state (Eqn. [l8| in H^y,r) with 7 ={0.33, 0.66, 
0.99}. The experimental results (Fig. |3(b)| ) shows excel- 
lent entanglement preservation and good agreement with 
the theoretical simulation. The experimental deviation 
of concurrence (Eqn. 19) is less than 2%. 
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CONCLUSION 



tanglement preservation mechanism suggested by Hou et 
al. 



We have performed Fidelity Profile Optimization for 
the Hamiltonian - DM interaction in the presence of 
Heisenberg XY interaction. The optimized UOD can be 
used for all relative strengths (7) of the interactions and 
length is invariant under 7 or evolution time. Using these 
decompositions, we have experimentally verified the en- 
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